Manifestations of a spatial variation of fundamental constants on atomic clocks, Oklo, 

meteorites, and cosmological phenomena 
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The remarkable detection of a spatial variation in the fine-structure constant a from quasar 
absorption systems [J must be independently confirmed by complementary searches. In this letter, 
we discuss how terrestrial measurements of time-variation of the fundamental constants in the 
laboratory, meteorite data, and analysis of the Oklo nuclear reactor can be used to corroborate 
the spatial variation seen by astronomers. Furthermore, we show that spatial variation of the 
fundamental constants may be observable as spatial anisotropy in the cosmic microwave background, 
the accelerated expansion (dark energy), and large-scale structure of the Universe. 



INTRODUCTION 

A very large study of quasar absorption systems has 
recently detected a spatial variation in the fine-structure 
constant, a — e 2 /Tic The general idea is to compare 
the wavelengths of atomic spectra measured in the lab- 
oratory with those seen in absorption systems at high 
redshift. Any variation in the value of a results in well- 
understood discrepancies between the two spectra. Sys- 
tematics are controlled to a very high degree by using 
many atoms and ions 0, Q . 

Previous studies of quasar absorption spectra had sug- 
gested that the fine-structure constant may have been 
smaller in the past [4j-l7(. However, these studies all used 
spectra taken at the Keck telescope in Hawaii, at a lati- 
tude of 20° N. Similar studies using the Very Large Tele- 
scope in Chile (latitude 25° S) at first showed a stringent 
null constraint [8]. More careful analysis of the same 
sample suggested that the errors should be enlarged by 
a factor of six 0, [13], and that a much larger, dedicated 
VLT survey be performed. 

The most recent work, which makes use of both Keck 
and VLT data, shows a highly significant (~ 4er) spatial 
gradient in the value of a [l|. That is, a was larger in 
the past in one direction and smaller in the past in the 
opposite direction. This dipole, which we will refer to 
as the "Australian dipole", has a declination of around 
—60°. This explains why the previous studies suggested 
a time- varying a that was smaller in the past: they only 
used data from the Keck telescope, which sees mainly in 
the Northern Hemisphere. The new results are entirely 
consistent with the previous ones. 

Discovery of a spatial variation in the fundamental con- 
stants of nature has massive implications for the "fine- 
tuning" problem. This is the question of why the con- 
stants of nature seem to be finely-tuned for life to exist. 
While the anthropic principle can be invoked to explain 
such tuning of the Universe, there remains the question 
of how it could come about. With the detection of spa- 
tial variation of constants we begin to have a natural 
explanation for fine-tuning: with many possibilities for 



combinations of constants all occurring within the (pos- 
sibly infinite in extent) Universe, we simply appear in the 
part of the Universe that is consistent with our existence. 

Extraordinary claims require extraordinary evidence, 
and the detection of a spatial gradient in a is no excep- 
tion. In this letter we discuss how limits from atomic 
clock measurements, meteorite data, and the Oklo nu- 
clear reactor may be interpreted in the light of the evi- 
dence for spatial variation. We show how these indepen- 
dent methods could be used to confirm or contradict the 
dipole model of a- variation that the quasar data suggest. 
We also briefly discuss the possibility that a spatial vari- 
ation may be observed via asymmetries in the expansion 
of the Universe, the cosmic microwave background, and 
the large-scale structure. 

MODEL 

The "Australian dipole" of a- variation found by [H is 



Sol 

— = (1.10 ±0.25) x lO~ 6 rcos-0 
a 



(1) 



where 5a/ ao = (a(r) — ao)/ao is the relative variation of 
a at a particular place r in the Universe (relative to Earth 
at r = 0). The function r cos?/> describes the geometry of 
the spatial variation: ip is the angle between the direc- 
tion of the measurement and the axis of the Australian 
dipole, (17.4(0.6) h, —58 (6)°) in equatorial coordinates. 
The distance function is the light-travel distance r = ct 
measured in giga-lightyears. This is model dependent for 
large redshifts: we use the standard Aqdm cosmology 
parametrized by WMAP5 [ll| to determine the light- 
travel time t. It is assumed here that Sa /ao = at zero 
redshift, which is supported by the data, however this 
assumption should be tested using the same absorption 
methods as are used at high redshift (e.g. by using ab- 
sorbers within our own galaxy). 

Our goal is to provide a simple interpretation of ter- 
restrial measurements of variation of constants in terms 
of the spatial gradient in values of a. A minimal hy- 
pothesis is to expect all fundamental constants to vary 
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along the direction specified by the Australian dipole 
of Ref. [l|. There are some good theoretical justifica- 
tions for this postulation. For instance, the constants 
may vary because they are coupled to a (dimensionless) 
scalar held $ which varies over space-time. For exam- 
ple the quintessence held $ /c 2 or a dimensionless dilaton 
field. In this case the axis of the dipole is the direction of 
its gradient V<fr, and a fundamental constant X is cou- 
pled to its variation via 

f) X 

-rr- = kxd$ (2) 

where kx is a dimensionless coupling coefficient. Our 
dipole model now requires <5$(r) ~ rcosip but all con- 
stants will vary in the same direction (i.e. along the 
Australian dipole). 

In this letter we will deal with the constants a = e 2 /he, 
the electron-to-proton mass ratio fi = m e /m p , and the 
dimensionless mass ratio X q = m q / Aqcd where m q is the 
light-current quark mass and Aqcd is the position of the 
Landau pole in the logarithm of the running strong cou- 
pling constant, a s ~ l/ln(AQCD^*/^c). In the Standard 
Model the electron and quark masses are proportional to 
the vacuum expectation of the Higgs field, while the pro- 
ton mass m p is proportional to Aqcd (if we neglect the 
~ 10% contribution of the quark masses). Relative vari- 
ation of X q is then approximately equal to the relative 
variation of m q /m p and fj, — m e /m p . We can relate the 
relative variation of different constants by equations like 

kfj, = k a (3) 

where the can be determined from observations and 
compared with theories of spatial variation. 



where tp is the angle between the motion of the Sun and 
the dipole. The best fit value for the Australian dipole 
gives cosV> ~ 0.07, but this has an uncertainty ~ 0.1. 
This signal will be modulated by the annual motion of 
the Earth around the Sun; with the angle between the 
ecliptic plane and the Australian dipole taken to be 35°, 
this modulation is 

— = 1.4 x 10- 20 coswt (5) 
a 

where uj refers to the angular frequency of the yearly 
orbit. 

Eq. (U) gives a useful benchmark for comparing labora- 
tory measurements to the recent spatial a- variation data. 
The current best limit on the rate of a- variation in labo- 
ratory measurements is a/a = (— 1.6±2.3) x 10~ 17 yr _1 , 
obtained by comparison of Hg + and Al + clocks over the 
course of a year [12| using the proposal and calculations 
of [1, [H . Therefore this limit will have to be improved by 
two orders of magnitude to compete with the astronom- 
ical spatial-variation data. Fortunately, atomic clocks 
are improving rapidly, and there are several new schemes 
that could allow measurement at this level of precision, 
including the 229 Th nuclear clock proposal 11311 that has 
a hugely enhanced sensitivity to a- variation [14|. 

To avoid misunderstanding, we note that the astro- 
nomical observations do not exclude temporal variation 
of a below the rate 10 -16 per year. That is, when con- 
sidering temporal variations, the laboratory observations 
are already competitive at their present level of accuracy. 



OKLO 



ATOMIC CLOCKS 

Laboratory-based atomic clocks provide measurements 
of the time-variation of a in the Earth frame. Since the 
Earth is moving with respect to the rest frame of the 
cosmic microwave background (CMB), we may expect a 
to vary as the Earth moves along the Australian dipole. 
The assumption here is that the dipole is stationary in 
the CMB rest-frame, which is to be expected if the varia- 
tion is caused by a co-moving scalar field, for example. To 
start, we neglect the annual motion of the Earth around 
the Sun, which is averaged out over time. The velocity 
of the Sun in the CMB rest-frame is known to high accu- 
racy from the CMB di pole itself, and is 369 kms -1 in the 
direction (168°, —7°) [11]. This is almost perpendicular 
to the direction of the Australian dipole, and is therefore 
sensitive to the exact angle, which has some uncertainty. 
We can expect the yearly variation in laboratory mea- 
surements to be 

d/a|iab = 1-35 x 10~ 18 cos ip yr" 1 (4) 



We can get an estimate for the kind of variation ex- 
pected at the Oklo natural nuclear reactor by seeing how 
far our local galaxy cluster has moved relative to the rest 
frame of the cosmic microwave background. Currently 
the center of our galaxy is moving at 552 kms -1 in the 
direction (266.5°, 29.1°) in galactic coordinates or (158°, 
—24°) in equatorial coordinates [l^. Extrapolating this 
velocity over the ~ 1.8 billion years since the Oklo reac- 
tor was running gives a total movement of 3.3 x 10 6 light 
years. Of course, this does not account for any partic- 
ular accelerations of the Milky Way galaxy with respect 
to the local group. Indeed, the local group as a whole is 
moving at around 630 kms -1 in roughly the same direc- 
tion: this gives an indication of the expected deviation 
from constant velocity. 

The galactic center moves at an angle cost/' ~ 0.24 to 
the direction of the Australian dipole, so it has moved a 
distance of 0.78 x 10 6 light years along the dipole direc- 
tion since the time of Oklo. Therefore, if the quasar 
results are correct, we expect that a will have been 
smaller on Earth in the past. To achieve accuracy from 
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Oklo comparable to the quasar data will require deter- 
mination of a-variation at the time of Oklo at the level 
Sa/a ~ 1.1 x 1(T 6 x 0.78 x 1(T 3 - 10~ 9 . The most sensi- 
tive Oklo limits on variation are based on measurement of 
the position of a very low energy resonance (E r =0.1 eV) 
in neutron capture by the 149 Sm nucleus. Using the re- 
sult [3 



5E r _ £Xg 

MeV ~~ X„ 



8a 
a 



(6) 



we expect that E r would have been different at the time 
the Oklo reactor was running by 



5E; 



r Oklo 



1 meV 



R% 10 meV 



(7) 



where R q is defined by k q = R q k a as in Eq. The 
Sm resonance has an order of magnitude more sensitiv- 
ity to variation in the dimensionless light-quark mass X q 
than it does to a. The current best 2a limits from Oklo 
data are -12 < SE r < 26 meV [17] and -73 < SE r < 
62 meV 0. 



in Gyrs. As in the Oklo case, this movement is that of 
the galaxy as a whole. Then 



^ = 4.8xio- 10 x?; 



and integrating ([8]) with to — ti = 4.6 Gyr gives 



A = Aq (l - 1.1 x 10~ 9 i^ 



(9) 



(10) 



where Ao = X(to). Therefore the relative difference be- 
tween the decay rates measured from meteorites and that 
measured in the lab is expected to be at the level 



A-Ap 
Ao 



= -2.2 x 10~ 5 



(11) 



Unfortunately, this is far below the current accuracy 
~ 10~ 2 , which comes from A = 1.666 (9) x 10~ n y r" 1 [23 



and Ao = 1.68(5) x 10~ n yr" 1 [2J] (see also [25j). Both 
of these measurements will need to improve significantly 
before the expected range of variation is reached. 



METEORITES 



COSMOLOGICAL IMPLICATIONS 



Meteorites can be used to determine average rates of 
decay of long-lived isotopes over the last 4.6 gigayears, 
which can be compared with laboratory rates to test for 
changes in fundamental constants [2(| ■ The energy 
difference in the /3-decay of 187 Re to 187 0s is very small, 
u) = 2.66 keV. It is the result of cancellation between the 
relatively large Coulomb and asymmetry energy differ- 
ences (which are ~ 20 MeV [19|, [2l|) of the two nuclei, 
and so it is very sensitive to possible variation of con- 
stants. 

We define our dimensionless observable by y = 
Afle/A[/: the ratio of the decay rate of the 187 Re iso- 
tope to that of a relatively insensitive isotope, such 
as Uranium. The relative variation of y is Sy/y = 
^Afle/A/fe — SXu/Xu- Since Xr c <C Xjj, the Uranium iso- 
tope serves to calibrate the time that the meteorite was 
formed. Therefore, taking Xjj as fixed, we can express the 
sensitivity of A^ e to variation of a as 8X/X = Sa/a, 

with k% ~ 2 x io 4 mm- 

The decay rate determined from the Re-Os 
isochrons 23] is an average value over the time 
since the meteorites were formed, 



A = 



to — t\ 



X(t) dt , 



(8) 



where t\ and to are the meteorite time and present 
time, respectively. As the solar system moves through 
the CMB rest-frame along the Australian dipole, we 
expect the fine-structure constant to vary as (JTJ) with 
rcosip — vt — 4.3 x 10~ 4 t where t and r are measured 



Variations of the fundamental constants will lead to 
changes in the masses and binding energies of elemen- 
tary particles, including leptons and baryons as well as 
dark matter particles. If one assumes cosmological con- 
servation laws apply in a flat Universe, then 



baryons 



9. 



dark matter 



radiation 



n 



dark energy 



= l 



and there must be an exchange between the light and 
dark matter energy densities and the "dark energy" . 
Therefore a cosmological anisotropy in the fundamental 
constants will result in an anisotropic energy exchange 
between these contributions to the total energy density 
of the Universe. 

The contribution of baryonic matter to gravitational 
processes will be approximately proportional to the di- 
mensionless constant Gm^/Tic ~ GAq CD /?ic. In the ab- 
sence of a confirmed theory for dark matter, we may 
write a similar expression for the dark matter contribu- 
tion, Grn^ )M /hc, and a cross-term GuipUiDyi / he, which 
depend on an unknown mass todm- The effect of varia- 
tion of these contributions to the energy density of the 
Universe may be related to varying G theories. 

We see that the observed a-variation is related to mass 
density and the cosmological "constant", and hence the 
accelerated expansion of the Universe. Therefore it may 
be possible to see a spatial anisotropy in the redshift- 
luminoscity relationships of Snla supernovae data. Ad- 
ditionally the temperature anisotropy of the cosmic mi- 
crowave background and large-scale structure formation 
may have observable spatial anisotropies. 
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CONCLUSION 

We have shown how laboratory experiments of time- 
variation of fundamental constants and limits from the 
Oklo nuclear reactor and meteorites can be compared to 
the spatial variations seen by astronomers. This interpre- 
tation is essential if one wishes to independently corrob- 
orate the spatial variation seen by astronomers. None 
of the current terrestrial limits contradict the reported 
observation of a spatial gradient in a. 

We have also tested other astrophysical systems where 
one could find spatial variation of fundamental constants: 
full results will be deferred to a later paper. Briefly, the 
existing Hi data from quasar absorption spectra shows 
hints that there may be a dipole in /i-variation with an 
axis corresponding to the a- variation of Ref. [l( . On the 
other hand the variation of x = a 2 fig p , inferred from 
21-cm data, has a best-fit dipole whose axis does not 
correspond to that of the Australian dipole, although in 
this case systematics heavily dominate. It is also pos- 
sible to infer spatial variation of fundamental constants 
during big bang nucleosynthesis from high-redshift mea- 
surements of primordial abundances. The existing deu- 
terium data does not support the dipole interpretation 
statistically, but the preferred axis is aligned with the 
Australian dipole. There is a strong impetus now to per- 
form measurements of relative primordial abundance at 
high redshifts of as many elements as possible in as many 
different spatial directions as possible. 

Finally we note that it may be possible to observe a 
spatial dipole in other cosmological systems. For exam- 
ple, a- variation may be seen in the CMB anisotropy if a 
high-enough sensitivity can be reached. Although the re- 
sults of [H (interpreted as strictly spatial variation) sug- 
gest that accuracy at the level 10 -6 will be required, if the 
hints from big bang nucleosynthesis turn out to be real 
then there is an additional redshift (time)-dependence 
that could increase the variation at the time of the CMB 
substantially. Another possibility is that since the ob- 
served a- variation is related to variation in the mass den- 
sity and cosmological constant, and hence the accelerated 
expansion of the Universe, it may be possible to see a 
dipole in the redshift-luminoscity relationships of Snla 
supernovae data. 

We thank G. t'Hooft, J. A. King, J. K. Webb, and M. 
T. Murphy for useful discussions. This work is partly 
supported by the Australian Research Council. 
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